Diagrammatic approximations for the 2d quantum antiferromagnet: exact projection 

of auxiliary fermions 
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We present diagrammatic approximations to the spin dynamics of the 2d Heisenberg antiferro- 
magnet for all temperatures, employing an auxiliary-fermion representation. The projection onto 
the physical subspace is effected by introducing an imaginary-valued chemical potential as proposed 
by Popov and Fedotov. The method requires that the fermion number at any lattice site is strictly 
conserved. We compare results obtained within a self-consistent approximation using two different 
auxiliary-particle projection schemes, (1) exact and (2) on average. Significant differences between 
the two are found at higher temperatures, whereas in the limit of zero temperature (approaching 
the magnetically ordered ground state) identical results emerge from (1) and (2), providing the 
qualitatively correct dynamical scaling behavior. An interpretation of these findings is given. We 
also present in some detail the derivation of the approximation, which goes far beyond mean-field 
theory and is formulated in terms of complex-valued spectral functions of auxiliary fermions. 

PACS numbers: 75.10.Jm, 75.40.Gb, 67.40.Db 
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I. INTRODUCTION AND REVIEW 

Auxiliary particles are a widely used tool in the theory 
of correlated electron systems. The principal difficulty in 
the treatment of these systems is the strong Coulomb re- 
pulsion U for two electrons on the same localized orbital, 
usually of d or / character. In effective model Hamilto- 
nians like the Heisenberg, t-J , or Kondo model the large 
U leads to a Gutzwiller projection onto the quantum- 
mechanical subspace, where none of the d- or /-orbitals 
may contain more than one electron at a time. 

In this paper we will focus on the spin-1/2 quantum 
Heisenberg antiferromagnet in two spatial dimensions 
(2D) on a square lattice^, 



(1) 



<ij> 



The sum covers all nearest-neighbor pairs < i,j > . The 
model (H|) may be obtained from the single-band Hubbard 
model for large U with nearest-neighbor hopping ampli- 
tude t (leading to J = 2t 2 /U) in the limit of a half-filled 
band with one electron per site 2 . It represents the sim- 
plest low-energy model for two-dimensional Mott insula- 
tors, in particular the CuO-planes in the undoped parent 
compounds for high-temperature superconductors^^. 

The restriction on states with no doubly occupied sites 
is reflected in the non-canonical commutation relations of 
spin operators, 



[S*,S V ] = iS z ^ c-number. 



(2) 



For an analytical approach to the model, Eq.(j2j poses a 
severe difficulty, since the standard Wick's theorem and 
many-body techniques cannot be applied 6 ' 7,8 ' 9 . A conve- 
nient way of circumventing this difficulty is to represent 
the spin operators in terms of canonical auxiliary-particle 
operators, of either fermionioA^ or bosonioii character. 
The cost of this concept is the extension of the Hilbert 



space into unphysical sectors. These unphysical states 
have to be removed by imposing a constraint. In this 
work we use auxiliary fermions, 



1 



2 ^ 



4>T a / ia , Qi /fa/fa = 1 , (3) 



, [i — x,y,z are the Pauli matrices, and a =|, j is 
the fermion-spin index. Here and in the following we let 
h = 1 . The representation ([3|) fulfills the commutation 
relations @ . The Fock space of the auxiliary fermions 
fi a is spanned by the states 

physical: | T) = , U) = /J|0) (4a) 
unphysical: |0) , |2) = /J/JlO) (4b) 



where |0) denotes the vacuum, fi a \Q) = . 

The projection onto the physical subspace can be 
performed in several ways 8 , and for impurity models 
like the Kondo model or the single-impurity Anderson 
model the projection of auxiliary particles is a standard 
tectmiqu o 12 ' 13 . (For the latter model, auxiliary bosons 
have to be added to the scheme^.) However, in lattice 
models like Eq.{T]) the constraint Qi = 1 has to be en- 
forced on each site i independently. This leads to the 
so-called excluded-volume problem 7 -^ and prohibits an 
infinite-order resummation of the perturbation series in 
J . In the limit of high spatial dimension the problem can 
be evaded within the (extended) dynamical mean-field 
theory. Here the infinite lattice system is approximately 
mapped onto a single sit o 16 ' 17 i 18 or a small cluster— cou- 
pled to a bath. 

An alternative path starts from a mean-field-like treat- 
ment of the constraint, where the auxiliary charge Qi = 1 
is fixed merely on the thermal average 



,,20 



Qi - (Qi) = (Qi) = J2(fUi*) = i 



(5) 



2 



This condition is introduced into the Hamiltonian {TJ, 
([5]) through a chemical potential pJ for the auxiliary 
fermions. Due to the particle-hole symmetry of ([1]) we 
have \J — . The approximation ([5]) is of great advan- 
tage, since now the perturbation theory in J starts from 
non-interacting fermions, and we can make use of the 
standard Feynman-diagram techniques. Eq.(J5|) is also the 
starting point for numerous mean-field theories of corre- 
lated electron systems. An improvement of the mean- 
field-like constraint ([5]) in a perturbative fashion has fre- 
quently been made by generalizing to a fluctuating 
Lagrange multiplier (see, e.g., Ref. @) ■ 

The Popov-Fedotov approach: The method proposed 
by Popov and Fedotov 21 enables one to enforce the 
auxiliary-particle constraint exactly within an analyti- 
cal calculation for the infinite syste m 22 i 23 i 24 i 25 . The ap- 
proach starts from a "grand-canonical ensemble" , 

i a 

(6) 

with a homogeneous, imaginary-valued chemical poten- 
tial 



f 7T 7T 

P J = t^k B T = i-T 



(7) 



H is the spin Hamiltonian 1[T]). written in terms of the 
auxiliary-fermion operators ([3]) . 

There are two main requirements for the Popov- 
Fedotov method to work: The first is the conservation 
of the auxiliary charge Qi on each lattice site, 



[Qi,H]=0 , i = 1,2,. 



,N T 



(8) 



where Nl denotes the number of sites (i.e., spins) in the 
system (JTJ . Since also [Qi,Qj] = , the eigenstates of 
H and H ppv in the Fock space of the fermions can be 
specified by some auxiliary-charge configuration 



cq = (<2l,<22, • 



) , Qi G {0, 1,2} 



(9) 



Physical states belong to the subspace with the configu- 
ration 



eg"" = (1, 1, .... 1) . 



(10) 



For a given configuration cq the Hamiltonian ([T]) has a 
specific set of eigenstates with quantum numbers uq and 
energies E{cq,tiq) , i.e., Schrodinger's equation reads 



h \cq, n Q ) = E(c Ql n Q )\cQ, n Q ) 



(11) 



Consider now the partition function Z ppv for the grand- 
canonical Hamiltonian ([S]), 

Z ppv = Tr^e-^H (12) 



Tr-^ denotes the trace in the enlarged Hilbert space of 
the auxiliary fermions, and j3 = l/fc^T = 1/T. With 
Eqs. pT|l and © it becomes 



Z ppv — 



E 



E 

=0 no 



• (c Q ,n Q \ e -^( c «'"«) e 1 ^ ■ ■ ■ e**«^ \c Q , n Q ) 

In addition to the Qi-conservation, the Popov-Fedotov 
method requires that the Hamiltonian and the opera- 
tors appearing in physical (i.e., observable) correlation 
functions destruct the unphysical states |0) , \2) . In the 
present case, Hamiltonian and correlation functions are 
composites of spin operators, and we have, using Eqs.((3|) 
and (Hbl). 



Sm = , |2) = 



(13) 



Consider an arbitrary site / with an unphysical auxiliary 
charge Qi ^ 1 . From Eq. (fT3")) it follows 

Qi^l: E{c Q ,n Q )\ Q[=0 = E(c Q ,n Q )\ Q _ 2 , (14) 

that is, the spin at site I seems to be removed from the 
Hamiltonian for all states \cq,tiq) with Qi = or 2 . 
Accordingly, the contribution from Qi — 0, 2 to Z ppv is 
proportional to 



Q,=0,2 



(1 + e i7r ) = 



(15) 



In that way, the unphysical contributions from all sites 
I = 1 , . . . , Nl cancel in the grand-canonical partition 
function, i.e., only the physical charge configuration 



c^ ys , Eg. ([TO]), survives in Z ppv 



Z ppv = {i) NL J2{n\e 



n) 



Here \n) and E n denote the eigenstates and -energies of 
the Hamiltonian (UJ in the physical subspace, 



\n) = |c^ ys ,n Q ) 



En 



E{^ vs ,n Q ) 



H\n)=E n \n) . 
Thus we end up with 

Z ppv = (i) NL Z , 



(16) 



(17) 



i.e., up to a constant prefactor, the (canonical, physical) 
partition function Z for the Heisenberg model (JTJ) is given 
by Z ppv . 

The above argument, originally presented by Popov 
and Fedotov in Ref. l2ll. is extended to Green's functions 
in Appendix [A]. It is found that any correlation func- 
tion of spin operators may be calculated from the grand- 
canonical Hamiltonian ([B]). In particular, the imaginary 
time-ordered spin susceptibility 



Xf(r-r') 



7 Tr[e-f )H T T {Snr)Sf(r')}[ 



(18) 



3 



can be obtained from 



Xlf (r - r') = (T T {S? (r)5f(r')}) ppv (19) 

with fJ,,p,£ {x, y, z} . The expectation value is calculated 
in the enlarged Fock space, 



(...)P pv = —Tr } \e-P HP 



(20) 



with Z ppv as dehned in Eq.([T2) above. The "modified 
Heisenberg" picture 6 - for spin operators reads 

S<t{r) = e TH S<?e- TH = e^Sf e"^ , (21) 

using again Eq. (|13p . Similarly, the local magnetization 
is given by 

(S?) = ^Trl e- pH S? } = (Sf ) ppv . (22) 

It should be emphasized that expectation values of 
wnphysical operators are meaningless within the Popov- 
Fedotov scheme: E.g., the auxiliary- fermion charge Qi 
introduced in Eq.j3]) docs not dcstruct the unphysical 
states (I4bl). and one has 



(23) 



For the l.h.s. we know that Qi = 1 in the physical Hilbcrt 
space. For the r.h.s., however, we obtain (Qi) ppv = (1 + 
i) . The calculation is given in Appendix [A"l. 

Average projection: For comparison, we also want to 
use the mean-field-like treatment of the auxiliary-fermion 
constraint mentioned below Eq. ([5]) . Since the (real- 
valued) chemical potential added to the Hamiltonian fl} , 
([3]) turns out to be zero, due to the particle-hole symme- 
try of ([T]), the calculation of a spin-correlation function 
or magnetization amounts to just enlarging the Hilbcrt 
space into the Fock space of the auxiliary fermions. The 
equivalents of the Eqs.lO, and (J23|) then read 



and 



whereas 



with 



(Qi) = {QiY 



(24) 



(25) 



(26) 



(■•■) av = ~n — Tr^[( 



-0H 



Tr^e-^1 



(27) 

The ~ sign in (|24|) , ([25]) stands for the error introduced by 
the uncontrolled fluctuations of the fermion occupation 
numbers Qi into unphysical states. These fluctuations 
are absent in the Popov-Fedotov scheme. 



In the following sections, results from the exact and the 
averaged constraint will be obtained within the same di- 
agrammatic approximations. The effect of the constraint 
on physical quantities like the dynamical structure factor 
and magnetic transition temperature is going to be stud- 
ied. It will turn out that at sufficiently low temperature 
the results for averaged and exactly treated constraint 
become equal. In addition we will show in some detail 
how a self-consistent approximation that goes far beyond 
mean-field theory, can be worked out within the Popov- 
Fedotov approach. 



II. EFFECT OF THE CONSTRAINT: SIMPLE 
APPROXIMATIONS 

In order to compare results for the Heisenberg model 
(TTJ from average projection, Eq.([5]), and exact projection 
using the Popov-Fedotov scheme, Eq.®, we consider a 
more general grand-canonical Hamiltonian of auxiliary 
fermions, H = H — \J J^i Qi ■ With the model Hamilto- 
nian ([T|) written according to ([3]) , it reads 



H 



(28) 



1,3 



Sums over spin indices a, a, (3, (3 are implied. The antifer- 
romagnetic coupling Jjj is nonzero and equal to J > , if 
i,j are nearest neighbors. Depending on the projection 
method, the chemical potential yJ takes the value 



average projection: [i 



f 



exact projection: 



/ 



, 

i-T 
2 



(29) 



In the case of exact projection, the Hamiltonian Ip8j) is 
no longer Hermitian. Nevertheless, physical quantities 
like the dynamical structure factor for spin excitations 
or the magnetization will come out real-valued. 

Eq. I|28p represents a system of canonical fermions with 
a two-particle interaction ~ J and may therefore be 
treated using standard Feynman-diagram techniques 6 . 
The bare fermion-Green's function, written as a matrix 
in spin space, reads 



1 



■ + uf 



(30) 



where u> = (2n + l)irT , n € Z is a fermionic Matsubara 
frequency. For the case of exact projection, \rf = iirT/2 , 
there is some common practic o 21 i 23 ' 24 to absorb \J into 
iui and to re-define u> accordingly. Here we do not follow 
this line, but keep lo as introduced above (i.e., fermionic). 

Free spins (J — 0): The simplest case is given by 
setting J — . Using Eq. J3|) , the susceptibility (|T8|) , to 
be calculated from Eq. lfT^) or (f2"4"|) , is then given by a 



4 



simple bubble, 



For average projection, with \J = , we find 



with, using Eq. l|3"Uj) . 

1 



X 



2T 



v = 2mrT is a bosonic Matsubara frequency, Tr CT denotes 
a trace in spin space, and f(x) — l/(e x ^ T +l) is the Fermi 
function. Depending on the constraint method, the result 
is 



average: /(0) = - x° 

exact: /(±4 T ) = ^ X° = 

2 1 ± « 



1 

8T ' 
1 

4T ' 



(31) 



Two observations are in order: The imaginary chemical 
potential cancels out in the physical quantity x° > an d the 
result with and without use of the exact auxiliary-particle 
projection is qualitatively the same (Curie law). With 
merely average projection in effect ({i* = 0), the spin 
moment S(S + 1) extracted from the Curie law x (T) <x 
S(S + 1)/T is reduced by a factor of 2 compared to the 
exact result. This is due to fluctuations of the fermion 
charge Qi . 

Mean-field approximation: The simplest approach to 
the interacting system J > is the Hartree approxima- 
tion, i.e., magnetic mean- field theory. This approxima- 
tion does locally conserve the auxiliary-fermion charge 
Qi . Dyson's equation for the auxiliary fermion reads 



Gi(iw) = [iu> + /jS — T,i(ito) 



(32) 



and in Hartree approximation the self energy is indepen- 
dent of iui and given by 



S * = Ji in a Tr<T t o a G J ( iuJ ^ e 



(33) 



(Si) 



Here the mean magnetization (Sj) on the site j has been 
identified. For a square lattice with coordination number 
z = 4 and only nearest-neighbor interaction J > we 
assume a Neel state on the two sublattices A, B , 

(8 A ) = -(8 B ) = -(S' B )e M . 
The fermion Green's function for any site on A becomes 



Ga(«w) = [iu + + o- z h] 



h= Z -±{S* A ), 



which leads to the self-consistent equation 

(S%} = kfi-h-n')-f(h-nt)} 



x h 

average projection: h = — J tanh (^jl) • (35) 

Within the Popov-Fedotov scheme, using \J = i^T , one 
has 

z h 

exact projection: h = — J tanh (— ) , (36) 

where the following expression for the Fermi function has 
been utilized, 

for a real-valued x . In the physical observable the 
imaginary part again cancels. Both projection schemes 
lead to the same self-consistent equation for the effective 
magnetic (Weiss) field h , except for a factor of 2 in the 
temperature. Accordingly the equations result in differ- 
ent Neel temperatures, 



average projection: T/v = —J 
exact projection: T/v 



However, at zero temperature both projection methods 
lead to the same result, 

average and exact projection: Hm (S^) = — ■ (38) 

That is, the unphysical reduction of the spin moment 
observed for the free spin, is completely restored in the 
magnetically ordered ground state. Apparently, the un- 
physical charge fluctuations are suppressed at T — » . 



III. EFFECT OF THE CONSTRAINT: 
SELF-CONSISTENT THEORY 

The purpose of this section is to demonstrate the ap- 
plication of the Popov-Fedotov scheme within a self- 
consistent approximation that goes far beyond mean 
field. To our knowledge, the Popov-Fedotov approach 
has at present been applied in mean-ficld-likc calcula- 
tions with perturbative correction o 21 ' 23 ! 24 ' 25 , but a self- 
consistent re-summation of the diagram series has not 
been attempted. 

When choosing an approximation scheme, it has to be 
kept in mind that the automatic cancellation of unphys- 
ical states requires the auxiliary-fermion charge Qi to be 
conserved (this has been discussed in Section U above) . 
In particular, the local gauge symmetry of the Hami- 
tonian (|28p under fi a 



(34) 



• fia must not be broken. 
Accordingly, approximations leading to finite expecta- 
tion values like (f} a fj a ) ^ or (/if/i|) cannot be 
used, while it is safe to consider so-called ^-derivable 
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FIG. 1: The self-consistent approximation (SCA) discussed 
in Section HU . Eqs. (|4Up . Continuous lines denote the fermion 
Green's function G , Eq. (|39[) , the thin dashed line represents 
J , and dots are Pauli matrices x 1/2 . 



approximations 26,27 . Spontaneous breaking of physical 
symmetries (e.g., spin rotation, lattice translation) may 
be included, since the respective order parameters are 
gauge invariant. As a consequence of the local gauge 
symmetry, the fermion Green's function is always local, 



Gi 3 (iuj) = 5ij Gj(wj) 



(39) 



Here we focus on the physics of the Heisenberg model 
([T|) in strictly two spatial dimensions at finite tempera- 
ture T > . This system has been studied extensively 
in the past using a variety of numerical and analytical 
methods, in particular in view of experiments on cuprate 
superconductors in the undoped limited. The Hartree ap- 
proximation discussed in Section |TT] above is, of course, 
not appropriate for the 2D system, although the Neel- 
ordered ground state at T = appears to be qualitatively 
correct^. At finite T the theorem of Mermin and Wag- 
ner requires the magnetization (Si) — to vanish. There- 
fore we seek an approximation, where the susceptibility 
X is self-consistently coupled back onto itself. Such an 
approximation has originally been proposed for the Hub- 
bard model^, and is commonly referred to as FLEX . For 
the Hamiltonian (|28|) with spin-spin interaction it takes 
the form shown in Fig.[T]. The fermion self-energy shown 
in the figure reads 



In the paramagnetic phase with lattice-translational sym- 
metry we have 

Gi(iuj) = cr° G{iuo) , %{iuj) = <j° E(iu) , 

and therefore 

The self-consistent equations (|40l) with Eq. (|3"2")l now turn 
into 

n(ii/) = ^2 G(iu + iu) G{kS) , (41a) 



x(q, iv) 



1 + J(q)II(ii/) ' 



(41b) 



H = TrE j2 ^)x^ w ) ' (41c 



D 



N L 
3T 



Y^D{iv)G{iu + iv) , (41d) 



G(iu>) — [ iu> + yJ — H(iu) ] 



(41e) 



The bare interaction in wave-vector space reads, for a 
square lattice in 2D with nearest-neighbor distance a = 
1, 

J(q) = 4J 7 (q) , 7 (q) = -[cosfe) + cos{q y ) } . (41f) 



%(iu) = ~TY Y Df{iv) a» Gi(iu + iv) a* , (40a) functional * close analogy^ 



Note that the bare J in Eq. (|40bj) does not contribute to 
the local D{iv) = Da(iu) , Eq. (|41c[) , since Ju = . 

It has been emphasized above, that the fermion prop- 
agator and as a consequence the irreducible bubble are 
local, Gij = Sij G , Hij = Sij n . Nevertheless, the inter- 
esting measurable^! quantity in the Eqs. (l4"T]) is the sus- 
ceptibility x(q) ! which is wave-vector dependent through 
the bare interaction J(q) , and therefore may describe 
even long-range fluctuations. 

The SCA shown in Fig. [1] can be derived from a 

to the FLEX. Therefore, 



uj and v denote a fermionic and bosonic Matsubara fre- 
quency, respectively, a Pauli matrix, /x, ft — x,y,z. 
The renormalized spin-spin interaction on lattice sites 
i,j is given by 



Dff H = -j y + ^JetffWJk. 



(40b) 



l.k 



The susceptibility x represents the series of fermion-linc 
bubbles in Fig. Q], 



xf = (n[i + jn] 



' I.) 



(40c) 



with 



nr(w) = -T^-Tr^a" Gi(iu + iv) a" Gj(iw)] 



(40d) 



Eqs. (l41|l represent a conserving approximation and can 
be used with the averaged fermion constraint as well as 
the Popov-Fedotov approach. The case of average pro - 



jection, fj,* = 0, has been treated in detail in Ref. |29| . 
The magnetic correlation length £(T) and the dynamical 
structure factor, derived from the susceptibility (|41b|) , 
came out quite satisfactorily when compared to known 
results, indicating that the diagrams in Fig. [1] indeed 
capture the important physics of the 2D system at low 
T . In the following some of the results from Ref. 
will be re-calculated using the Popov-Fedotov method, 
i.e., ju' = i\T . It will turn out, that the imaginary- 
valued chemical potential requires the use of complex- 
valued spectral functions, leadin g to more involved equa- 
tions than those derived in Ref. |29| for \J = . The re- 
sults calculated with both methods differ, except in the 
limit of vanishing temperature, where average and exact 
projection become equal, as will be presented below. 
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Equations for exact projection: Before discussing the 
numerical solution of Eqs. (|4ip , we quote some important 
formal results derived in Appendix El. For a numerical 
implementation it is suitable to express all Green's func- 
tions through their respective spectral functions. On ac- 
count of the Hamitonian being non-Hermitian, the spec- 
tral function of the fermion Green's function G{iui) be- 
comes complex-valued. G is given by 

(G.r a M= f 1/T dre^{T T {f la (r)fL(0)}f , (42) 
Jo 

with the thermal expectation value and r-dependencc 
calculated in the enlarged fermion Fock-space with the 
Hamiltonian (f2"5)) . G has the following spectral repre- 
sentation, G^ a {iui) — S a ,aG(iu)) , 



G(iuj) 



de- 



G(e) 



oo iuj + n f ~e 



^=i-T. (43) 



The energy variable e is a real number. For simplicity a 
system invariant under lattice translations and spin ro- 
tations has been assumed. The spectral function G is 
complex- valued , 

G(e) = Pl (e) + ip 2 (e) . (44) 
In the Appendix [A] we also derive the sum rule 



deG(e) = 1 



(45) 



=> Jdepi{e) = l , J dep 2 (e) = , 

e is again real valued. For the special case of spin de- 
generacy considered here, the fermion spectral-function 
in addition obeys the "particle-hole" symmetry, 



in detail in Appendix [Bl. Eqs. (|4ip now turn into 

J\f(e) e 2 



U(u) = 



de n o (,47a) 

oo (i + en'H) 2 + ( e n"H) 2 



II" (w) = S°(uj)-S (~uj) 
p ~ den"(e) 



n e — u> 



S°(lo) = - I <\:G + i:)G (.: - ^ ] . 



(47b) 
(47c) 

(47d) 



for the real (physical) functions U, S° ,H" ,H' with the 
Kramers-Kroenig transform II' , and 



G+M = [l-/( w -t-T)]G( W ) , 
G-{u) = f{u-i^T)G{u) , 



G{w) 



E( W ) 



E(w) 



( W -EH) 2 + (ttSH) 2 ' 



(47e) 
(47f) 

(47g) 



4tt 



deU{s) [G + (uj-e)+ (47h) 



+ G-{uj + e)~ 



w — e 



(471) 



for the complex (unphysical) spectra G ± ,G, E, E with 
the Hilbert transform E . Note again that the energy 
arguments u), e are real-valued. As is also shown in the 
Appendix [B], the Eqs. lfTT]) can be somewhat simplified 
by utilizing the symmetry (|46p . The imaginary chemical 
potential i^T appears in Eq. (|47ep and (|47f]) . adding an 
imaginary part to the fermion spectra G + and G~ via 
Eq.©. 

The density-of-states that enters Eq. (|47a|l is defined as 



A/ '( £ ) = ] v -E (5 ( e - J (q)) 



(48) 



G(-e) = G(eY 



Pi(-e) = pi(e) , p 2 {-e) = -p 2 (e) 



(46) 



The relations (gSJ), CH}, ([46]) hold as well for the fermion 
self-energy E . 

For a numerical solution of the Eqs. (|4Tj) , we introduce 
structure factors for Tl(iv) and D(iv) according to 

S°(li) = [ 1 + ] II"(w) . n"(w) = ImH(a; + iO + ) , 



U(aj) = [l + g(oj)]ImD(uj + iO + ) , 

with analytic continuation to the real axis via II (if) — > 
n(w + i0+) and -> D(w + i0+) ,u£l. The Bose 

function is denoted by g(uj) = \/{e u / T — 1) . As shown 



and for the nearest-neighbor interaction (|41fp in two di- 
mensions it becomes 



(49) 



with the complete elliptic integral of the first kind, 
K{m) = dt [ (1 - t 2 )(l - rnt 2 ) ]-^ 2 . 

The physical output from the numerical iteration of 
Eqs.([47|). (|49|) is the structure factor C/(w) of the effec- 
tive interaction, which is essentially the local (on-site) 
spin-excitation spectrum, and the wave- vector dependent 
dynamical spin-structure factor 



S(q,w) = [l+ fl (u;)]Im 



n( w ) 



1 + J(q)n(w) 



(50) 



with n = n' + ill" . Furthermore, the magnetic corre- 
lation length £(T) is extracted from the static magnetic 
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susceptibility x(q, 0) : Eq. (|41bp . for q close to the Neel- 
ordering vector Q = (tt, it) , i.e., 



1200 



leads to 

x(q,o) = 



J(q)~-4J+(q-Q) 2 J 

ri'(o) 



1 



i + j( q )n'(o) j^ 2 + (q-Q) 2 



where 11(0) = II' (0) has been used, and the correlation 
length is identified as 



jn'(o) 



1/2 



w=[T=um) (51) 

For completeness, in Appendix [C] the self-consistent 
equations (|47[) are re- written in real- valued spectral func- 
tions. These Eqs. (IC2p can directly be compared to the 
Eqs.(Al) in Ref.l29l: Both sets of equations represent the 
same diagrammatic approximation shown in Fig. [T], the 
former derived within the Popov-Fedotov scheme (exact 
projection), the latter within average projection. 



IV. NUMERICAL RESULTS 

The numerical results presented below are obtained 
from an iterative solution of the self-consistent equations 
(HIJ) , using the identical procedure and parameters for ex- 
act projection = i^T , leading to Eqs. (|C2p in App. 
[C| as well as average projection = , leading to 
Eqs.(Al) in Ref.lH). The procedure also utilizes Eqs. (|49"l) 
and (|5Tj) . which apply to both projection schemes. The 
data shown in Ref. [29| for the case of average projection 
have not been re-used in the present study. 

Magnetic correlation length: The correlation length 
£(T) is shown in Fig. |2]. £(T) becomes larger than one 
lattice spacing for T < J and reaches values up to ~ 1200 
lattice spacings for the lowest temperature T = 0.048J 
considered in this work. The numerical data shown in 
Fig. [2] are well reproduced by 



£(T) = c(T) 6 exp(a^) 



(52) 



The parameters a, b, c are determined by plotting the 
data as ln(£) vs. J/T (see the inset of Fig. [2J and a nu- 
merical fit of the function ln(c) + b ln( J/T) + a J/T to 
the data. We find 



average projection: a 
b 



exact projection: a 
b 



0.296 
0.592 
0.393 

0.304 
0.418 
0.595 



(53) 



At low temperatures the effect of the exact auxiliary- 
particle projection is only marginal, as is already ap- 
parent from inspecting Fig. |2]. The "spin stiffness" 
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FIG. 2: Main figure: Magnetic correlation length £(T) 
from the numerical solution of the self-consistent approxima- 
tion Fig.(fl| for the 2D square lattice, in units of the lat- 
tice spacing. The continuous line with diamonds corresponds 
to Eqs. (|47p . where the Popov-Fedotov approach has been 
used. The dashed line with circles belongs to Eqs.(Al) in 
Ref.|2^, where the auxiliary-particle constraint has been ap- 
proximated by its thermal average. Inset: The correlation 
length, plotted as ln(£) vs. J/T . Shown is the temperature 
region 0.048 J < T < 0.16 J , where a fit to the data has been 
performed as described in the text. 



a = 2irp s /J in the exponent in Eq. (|52p is the same in 
both cases^, a ~ 0.30 , merely the power b ~ 0.5 of 
the algebraic prefactor is slightly modified in going from 
average to exact projection. 

Spin spectral-function and energy scale: The almost 
vanishing effect of the auxiliary-particle constraint at low 
temperature is also visible in the spectra: Fig. [3] shows 
the effective interaction U(ui) , Eq. (|47a|) . which is the 
structure factor of the D(iv) given in Eq. (|41cp . Since 
J(q) 2 in Eq. (|41cp depends only weakly on q, U(uo) is 
essentally the local magnetic structure factor or spin- 
excitation spectrum. Note that Eq. (|47a[) is the same in 
the average-projected case, Eq.(Alh) in Ref. ,29. The 
data for low temperature, shown in the bottom panel of 
Fig. [3], features a broad shoulder of width ~ J , which 
is reminiscent of the box-like density-of-states for spin 
waves in 2D . Around zero energy U(lo) displays a huge 
peak (see the inset of the figure), which contains the 
critical fluctuations at q ~ (it, ir) close to the antifer- 
romagnetic ordering wave- vector. The curves for average 
and exactly treated constraint are on top of each other; 
merely the amplitudes of the peaks at ui = differ by 
a factor of 0(1) . At high temperature T = 2 J , on the 
contrary, the curves for average and exact constraint are 
well separated, in particular the total spectral weight is 
smaller for average projection. This can be seen from the 
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average constraint 

exact constraint 



T=2.0J 




2 3 
Energy co/J 

FIG. 3: Top: The structure factor U(u) , Eq. (|47a|) . of the 
renormalized local interaction D(iv) , Eq, (141c[l . at high tem- 
perature T = 2 J . Dashed and continuous lines as in Fig. [2]. 
Bottom, main figure: U (uj) at low temperature T = 0.05 J , 
corresponding to a correlation length £ ~ 900 . Bottom, in- 
set: Same data as in the main figure, but zoomed around 
the energy uj = . For clarity, the peaks for average con- 
straint (dashed) and exact constraint (continuous) are shifted 
from their original position uj = by 2e-05 J and -2e-05 J , 
respectively. 



data in the top panel of Fig. [3J. The reduction of spectral 
weight in U (uj) is related to an unphysical lack of local 
spin moment, which occurs if the constraint is not taken 
exactly. This will be discussed further below. 

The difference of peak amplitudes visible in the inset 
of Fig. [3J bottom can be traced back to the slightly dif- 
ferent correlation length £ (T) , compare Eq. {55J . The 
influence of the absolute value of £ vanishes if the scal- 
ing behaviour of U (uj) is considered: We start from the 
dynamical scaling hypothesis 33 , 



S(q,u) 



1 

UJQ 



S st (Q)<p(kOHH,uj/to ) 



(54) 




Energy <B/Cfl 

FIG. 4: Main figure: The scaling function Eq.j58} of the 
renormalized spin-spin interaction, using the energy scale 
uj (T) from Eq.JSTI and T(0) = 1. The curve shown here 
is the same for all temperatures 0.048 < T/J < 0.13, cor- 
responding to 1000 > £(T) > 10, and independent of the 
auxiliary-particle constraint being treated exactly or on aver- 
age. Inset: Detailed view of the region near u> = . 



Here S(q,u>) is the dynamical structure factor Eg . (f50|) . 
S st (Q) denotes the static (equal-time) correlation func- 
tion 



i r°° 

5 s * (q) = <S*S* q ) = - / d^5(q, W ) , 

ft J — oo 



(55) 



taken at the AF ordering wave vector Q = (tt, it) . tp(x) , 
&(x,y) are scaling functions, k = q — Q , and ui is the 
energy scale for critical fluctuations. At small energies 
ui <C J we expect U(ui) oc J d 2 q S(q, uj) , and from an in- 
tergration of Eq. ([54]) over wave- vector space there follows 
the scaling property 



U(uj) 



1 



S st (Q)T(Lu/uj ) . 



(56) 



T(y) is the (a priori unknown) scaling function for the 
local effective spin-spin interaction U (uj) . According to 
Eq. (f54"|) the energy scale can be extracted from the nu- 
merical data, up to a constant prefactor, using 



UJ 



S st (Q) 



S(Q,u) 



(57) 



j=0 



We obtained the energy scale in the temperature range 
0.048 < T/J < 0.13, which corresponds to correlation 
lengths 1000 > £ > 10 , using Eqs.([57|), J55J , and JSDJ) . 
The scaling function is then determined for each temper- 
ature from Eq. (l56|) . 



T(uj/ujo) 



U(0) 



(58) 
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FIG. 5: The spectral function G(u>) of the auxiliary fermion, 
introduced in Eqs. (|42|) . (|43|l . For exact auxilary-fermion pro- 
jection it is complex valued and given by Eq.( |47g| ), for average 
projection Eq.(Ale) from Ref. l29l holds, with G(u>) being real. 
Top: high temperature, Bottom: low temperature. Param- 
eters as in Fig. [3]. 



All curves T(w /u>q) , whether calculated with average or 
with exact auxiliary-particle constraint, agree to within 
numerical accuracy. The scaling function is shown in Fig. 
H]. Deviations from scaling behaviour become visible only 
for higher energies co/too > 20 . In particular, going from 
exact to average projection has no effect on the scaling 
behavior. 

A slight difference in the low-temperature properties 
of the two auxiliary-particle methods shows up in the 
energy scale itself: From a linear regression on u)q(T) x 
£ 2 (T) , obtained from Eq.([57| and Eq.([5"T]). we find 



average: wo£ 
exact: w £ 2 



(0.234 ±0.002) J 
(0.234 ±0.002) J 



0.089 T . 



(59) 



Nevertheless, the temperature behavior of the energy 



Temperature T/I 



FIG. 6: The local spin moment Sfo C as defined in Eq. (|61|l . 
multiplied by 4 . The limiting values for T J given in 
Eq. (|63[l are indicated as dotted lines. The exact value, given 
by the sum rule lf62]). is 45f* c = 1 . 



scale essentially is the same, 

wo(T) =0.23- J(£)- 2 , (60) 

which corresponds to a dynamical critical exponent z — 
2. 

Fermion propagator: Fig. [5] display the spectral func- 
tion of the auxiliary-fermion propagator for the two pro- 
jection methods. As has been discussed below Eq. (|4"2")l . 
the spectrum G(u>) of the fermion is complex valued, if 
the auxiliary-particle constraint is enforced exactly via 
the imaginary- valued chemical potential = i\T , If 
the constraint is treated on the average using [i? = , 
G{u>) remains real. In Fig.[5]the spectrum G{lj) is shown 
for high (top panel) and low temperature (bottom). At 
high T the spectra for exact and average constraint differ 
significantly, in particular the G(u>) from exact projection 
has a considerable imaginary part. At low temperature, 
on the other hand, the imaginary part is quite small, 
while the real part becomes equal to the spectrum of the 
average projection. In accordance with the results for 
spin-structure factor and correlation length, the projec- 
tion onto the physical part of the fermion- Hilbcrt space 
has almost no effect at sufficiently low temperature. 

Local spin moment: Another interesting quantity for 
studying the influence of the auxiliary-fermion constraint 
is the local spin moment, given by the local equal-time 
correlation function at an arbitrary site i , 



N, 



q 



(61) 



The static structure factor 5(q) has been defined in 
Eg. ([55)1 . For a spin- 1/2 system in the paramagnetic 
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phase Sf j* should fulfill the sum rule 



QSt 
D loc 



^=4 



(62) 



At very high temperature T J , interaction effects can 
be ignored, and Sf* c is given by the simple bubble II 
shown in Fig. [I], calculated with free auxiliary fermions. 
With Eq. (|41ap we have 



T> J 



QSt 

°loc 



iv 

^(/j-f/it) (filfil) \j=0 



For both projection methods, the expectation values are 
to be taken in the generalized grand-canonical ensemble 
([2"5| . (f2"5|) . For the case of average projection we find 

(fiafw) = /(0) = V 2 f° r an Y s P in direction <r = ±1 , 
with exact projection we have (fj a fi a ) = f{—i§T) = 
|(1 + i) . / denotes the Fermi function. That is, for 
T> J, 



average constraint: 4 



exact constraint: 



= 2 ||(1 = 1 



(63) 



4 <? s< 



If the projection onto the physical Hilbert space is per- 
formed exactly, the sum rule (l62|) is correctly reproduced. 
With average projection, on the other hand, it is signif- 
icantly violated. This is due to thermal charge fluctua- 
tions into unphysical, spinless states, which reduce the 
spin moment. 

When temperature is lowered, we find a partly unex- 
pected result: Fig. [6] shows 4 5^ c from the numerical 
solution as function of T . At high temperature the free- 
spin result is approached, whereas at low temperature the 
average and exact projection methods lead to the same 
value for Sf^ c . This observation fits into the line of results 
obtained so far: at T — > it does not matter whether the 
auxiliary-particle constraint is treated exactly or on the 
thermal average. 

However, the local spin moment at T — > , 
0.75 is too small. This is not due to an ill-treated con- 
straint, but an artifact of the approximation to the in- 
teracting system. The local moment (l61j) measures the 
total spectral weight of spin excitations, averaged over 
the Brillouin zone. The self-consistent approximation we 
use, see Fig. Q], apparently lacks some weight in the spin- 
excitation spectrum. In Ref. |2^ we studied an approxi- 
mation with a somewhat reduced self consistency (called 
"MSCA"), which delivered a better result for Sf* c at low 
T, iSi* c ~ 0.85. Moreover, the wave-vector depen- 
dence of S(q,ui) came out better. However, the MSCA 
cannot straight-forwardly be extended to the Popov- 
Fedotov scheme, since we don't know a ^-functional for 
that approximation to guarantee the conservation of the 
auxiliary-charge Qi , which is a necessary condition for 
the Popov-Fedotov method (see Section |l] above) . 

The limit T — > : At first sight it seems trivial that 



almost no effect at T — > : ^ enters the self-consistent 
equations through the Fermi function f(uj — i^T) in 
Eqs. ([47ej) and (j47fj) . Assuming T < ui , Eq.([37]) yields 



T«^: f(u>-i~T) = Q{-uj) + 0{e- 2 ^ T ) , (64a) 

which matches the case yJ — (average projection) at 
T = . On the other hand, the energy scale ojq of spin 
excitations, Eq. (|6"0"|) , is exponentially small compared to 
T , therefore the opposite limit should apply, 

T»c: /( w -i|T) = i(l+i) + 0(^) , (64b) 

adding a significant imaginary part to the fermion 
spectral-function G(io) . It requires a solution of the set of 
Eqs. (|47|) . however, to reveal that G{ui) has no features^ 
at w < wo (see Fig. [5]). The fermion spectrum is gov- 
erned by its bandwidth ~ J , and therefore the crossover 
from high-temperatures corresponding to Eq. (|64b[) to low 
temperatures, where Eq. (|64aP becomes valid, happens at 
T ~ J . A more physical interpretation is provided in the 
following Sections [V] and I VII . 



V. MEASURING THE CONSTRAINT 

In order to understand the weak influence of the 
fermion constraint at low temperature, it is useful to cal- 
culate the auxiliary-charge fluctuations in average pro- 
jection. Starting from Eq. (j2"Tj) , we have to calculate 



(AQi AQ,) av = (Q i Q J } Kr - (Q 4 ) av (Q 3 ) a 



(65) 



Since all charges Qi are conserved, [Qi , H] = , 
[ Qi , Qj ] = , the correlation function (|6"5j) may be ob- 
tained from the charge propagator in Matsubara space 
as 

Xg(i") = ^ T dr ^ iVT {%-{AQi{r) AQj(0)}) av 
Jo 

= i S v>0 (AQi AQj) av . (66) 

X^j(iv) is conveniently calculated with the Feynman- 
diagram rules introduced in Sect. [II], using the Hamil- 
tonain (|28[) with ^ = and the bare fermion Green's 
function (H . 

Free spins: In Section |TT] we first discussed the limit 
J = . In that case, is given by a simple bubble of 
bare fermion Green's functions, 

XtfW = -T£Tr CT [G°.^ + HGjM] 

i 



2T 



that is, 



the imaginary-valued chemical potential /i 



/ 



if T has 



J = 0: (AQiAQj) 



, av =<$■•- 



(67) 
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As expected, the auxiliary-charge fluctuations are finite. 
Note that X%{^ v ) is local (~ Sij) and static, (~ 5 V fl) in 
accordance with Eq. ([6"6")) , i.e., the conservation of the Qi . 

Mean-field theory: The second example presented in 
Section |TT] is the Hartree approximation, x® is again 
given by the simple bubble^, with £?y replaced by 

G y -M = 6 h3 [iu; + {-l) l a z h}- 1 , 

with the Weiss field h given in Eq.([33|) . For T <C J we 
obtain 

Hartree: (AQ t AQj) av = 8^2 exp(-/i/T) , (68) 

with h = ^ + <D(er J l T ) . That is, in the magneti- 
cally ordered phase the unphysical charge fluctuations are 
strongly suppressed, since the fermions develop a charge 
gap similar to a spin-density-wave state. 

Self-consistent theory: The approximation discussed 
in Sections IIIII and IIVI requires a more elaborate calcu- 
lation. For the case of average constraint, ^ — , the 
approximation given by Eqs.pT|) and Fig. [1] has been 
studied earlier in Ref . [23 ■ In Section IV C of that pa- 
per the conserving-approximation scheme has been ap- 
plied to the spin susceptibility, leading to a vertex func- 
tion in the fermion bubble, which is determined by a 
Bethe-Salpeter equation. The corresponding diagrams 
are shown in Fig. 6 of Ref.H^. For the auxiliary- fermion- 
charge susceptibility we want to calculate, the diagrams 
are exactly the same, except that the two spin vertices 
appearing in the bubble II in Fig. 6 are to be replaced by 
charge vertices er° = 1 . The response function (|66| then 
reads in wave- vector space 

X Q (q,0) = -T^Tr>°G( W )r(q, W )G( !W )] . 
For (AQ AQ) av merely the static limit lim x® (q, iv) is 

v— >0 

needed. With spin-rotation symmetry one has G = <j° G , 
and it turns out that T = cr°T (i.e., the charge and spin 
channels do not mix in the vertex function), leading to 

X Q (q, 0) = -2TJ2 GM 2 T(q, iu>) . (69) 

The vertex function is specified through the following 
Bethe-Salpeter equation, taken from the diagrams in Fig. 
6 of Ref. Hi, 

r(q,tw) = 1+ 

+ ^T^G(^i) 2 D{iu - iw 1 )r(q, iux)- 

iv k 

xTJ2 G{iuji) 2 [ G{iu3t + iv) + G(iu! - iv) ] T(q, 

10)1 

(70) 



The fermion Green's function G{iu>) and the local spin- 
spin interaction D(iv) have to be taken from the solution 
of the Eqs. (j4Tj) for \J — . The non-local spin interaction 
appearing in ([70)) reads 

D(q,zi/) = -J(q) + J(q) 2 x(q,^) , 

with x from Eq. (|41b|) . 

In Appendix iDlit is shown that the 2nd term in Eq. ([70[) 
actually becomes zero by symmetry arguments, i.e., the 
Bethe-Salpeter equation simplifies to 

r(q, iuj) = 1+ -T^2G(iui) 2 D(iu)-iui)T(q,iu>i) (71) 

Instead of solving the last equation numerically, we find 
it more instructive to aim at an approximate analytical 
solution. We employ the static approximation introduced 
in Ref. 29i, i.e., let D{iv) = D(0)S Uj0 in Eq.([7T]) as well 
as the fermion self-energy (|41d[) . The calculation closely 
follows Sections IV A and C of Ref. , leading to 

1 16 

1 - K/2)^ ' ^ = J ^ + °^ 

for temperatures T <C J ■ ojf is a typical 1/2 bandwidth 
of the continuous fermion spectrum, compare the bot- 
tom panel of Fig. [5]. Performing the Matsubara-sum in 
Eq. ([6"9"]) eventually leads to 

X Q (q,0) = — *«(T/J) , ^(t) = - + 0(t 2 ). 

Uf IT 

If the vertex function is ignored, T — > 1 , the result does 
not change significantly, § Q (t) -> || + 0{t 2 ) . Note that 
X® is independent of q, i.e., local. 

From Eq. ([6"6"]) we thus find the auxiliary-charge fluc- 
tuations of our self-consistent approximation in average 
projection, 

self cons.: (AQi AQj) av = 5 id — + 0(T 3 ) . (72) 

Since the fermion spectrum is gapless^l around u> = 
(see Fig. [5]), and the vertex function has very little effect, 
X® comes out Pauli-like. The explicit T factor in Eq. ([6"6"|) 
suppresses the charge fluctuations at low temperature. 
Compared to the magnetically ordered state described 
in mean-field (Hartree) theory, where the fermions de- 
velop a gap ~ J (see Eq. ([6"8")l ). the suppression of charge 
fluctuations is much weaker—. However, at T — > the 
unphysical fluctuations still vanish, and the average pro- 
jection becomes exact. 

VI. SUMMARY AND CONCLUSIONS 

Resummed perturbation theory is a powerful tool for 
calculating dynamical properties of strongly correlated 
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electron systems. In this paper we focused on the spin- 
1/2 antifcrromagnetic quantum Heisenberg model on 
the two dimensional square lattice. Summing infinite 
classes of contributions in perturbation theory is most 
economically done within a quantum-ficld-theoretic for- 
mulation employing canonical fields. We therefore use an 
auxiliary-particle representation of spin operators, which 
is a faithful representation in the physical sector of the 
Hilbert space. 

The use of auxiliary fermions requires a projection onto 
the physical part of the fermion-Fock space, where the 
fermion charge Qt equals one for each lattice site. 

While for a single lattice site the projection may be 
done exactly, e.g., by introducing an auxiliary-fermion 
energy A , which is sent to infinity at the end of the 
calculatio n 10 ' 14 , these standard methods cannot directly 
be generalized to effect the projection at each lattice site 
independently (this would require handling a large num- 
ber of independent limiting procedures, an impossible 
task in practice). 

For lattice systems the most simple approach to the 
projection is an approximative treatment, where a global 
chemical potential (Lagrange multiplier) [J is intro- 
duced, which is sufficient to fulfill the constraint on the 
thermal average, (Q,-) = 1 . 

However, Popov and Fedotov have proposed a rather 
unusual projection method, where a global imaginary- 
valued chemical potential /z* = i^T leads to an exact 
cancellation of unphysical states, therefore enforcing the 
operator constraint Q,- = 1 . Unfortunately, the Popov- 
Fedotov method may not straight-forwardly be general- 
ized to systems away from particle-hole symmetry 22 . 

In this paper we explored the usability of this con- 
cept by identifying the conditions to be satisfied by any, 
necessarily self-consistent, approximation scheme. Most 
important is the conservation of the fermion charge Qi 
by the model Hamiltonian, [Qi,H] — . If the approx- 
imation under consideration violates this conservation 
law, results become meaningless (see Sec. Hand App.|A| ■ 
Therefore, self-consistent approximations are most safely 
based on the conserving-approximation principle. Any 
hopping of auxiliary fermions, for example, is precluded 
by this requirement: The fermions are strictly local en- 
tities. The physically observable momentum dependence 
of spin correlation functions originates from the momen- 
tum dependence of the exchange interaction. 

Within the Popov-Fedotov approach the well known 
Feynman-skeleton-diagram expansion is applicable in 
conjunction with an exact projection of the auxilary 
particles onto the physical Hilbert space. We have 
shown in some detail how a self-consistent approxima- 
tion, which goes far beyond mean-field theory similar to 
the "fluctuation-exchange approximation" , can be for- 
mulated using complex-valued spectral functions of the 
(unphysical) renormalized fermion propagator. The re- 
sulting equations have been solved by numerical itera- 
tion. 

We applied the Popov-Fedotov method on several ap- 



proximation levels: the free spin, the Hartree approx- 
imation (magnetic mean-field theory), and the above- 
mentioned self-consistent approximation, using both av- 
erage projection (^ = 0) and exact projection (^ — 
i^T) . The results obtained for the latter approxima- 
tion show the expected suppression of the ordered state 
down to zero temperature, the exponential divergence of 
the spin correlation length, and a spin-structure factor 
consistent to the dynamical scaling hypothesis. 

A comparison of the results from average and exact 
projection reveals that there is a significant effect of the 
exact projection at higher temperatures. In the limit of 
low temperature, however, the deviation of the average- 
constraint results from the exact-constraint results be- 
come (numerically) indistinguishable, except for the case 
J = (free spins). 

In order to support this observation, we calculated the 
fluctuations (AQiAQj) of the auxiliary-fermion charge 
within the average-projection scheme. We find (by an- 
alytical calculation) lim^^o ( AQi AQ,, ) = 0, except 
for the case of free spins, where (AQi AQj) stays finite 
as T —> . That is, as long as the spin-spin interac- 
tion J is taken into account, the fermion-charge fluctu- 
ations into unphysical Hilbert-space states are quenched 
at T = . If temperature is increased from zero, we find 
that (AQiAQj) raises continuously with T . 

These at first sight surprising results find their expla- 
nation in the tendency towards antiferromagnetic order 
in the interacting system, which helps to suppress the 
fluctuations in the fermion-occupation number: Starting 
from the physical ("true") ground state, which features 
long-range magnetic order—, a fluctuation of the fermion 
charge Qi = 1 at some site I into an unphysical stated 
with Qi = or Qi = 2 is equivalent to removing the spin 
Si in the Hamiltonian (recall Eqs.Cd]) and (H])) . The 
lowest-lying state in this unphysical subspace thus lacks 
the binding energy of the spin at site I , which is of order 
,/ . Therefore, the ground state in the Fock space of arbi- 
trary fermion occupancy is the "true" ground state in the 
physical segment, and the lowest-lying unphysical state 
is separated from the ground state by a gap2& AEq ~ J . 

Consequently, at low temperatures T <C AEq , to a 
good approximation the exact projection may be omitted 
in favour of the technically somewhat simpler average- 
projection approach. At T = the approximate treat- 
ment of the constraint even becomes exact. Note that 
T <C AEq does not impose any restriction on the excita- 
tion energy ui , e.g., in the structure factor 5(q, lu) : Since 
the fermion charge is conserved locally, all excitations at 
any oj out of the ground state remain in the physical 
Hilbert space. 

The above argument is quite apparent for magneti- 
cally ordered systems. However, it should also apply to 
systems without magnetic order but strong correlations 
in the ground state. Examples are the various valence 
bond states discussed for, e.g., Heisenberg models with 
frustration^. In these systems the gap AEq to unphys- 
ical states is also expected to be ~ J . Somewhat dif- 
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ferent examples are systems with a ground state that is 
dominated by local Kondo singlets. Here the gap AEq 
is exponentially small in J , since the binding energy of 
a localized spin to the Fermi sea is given by the exp. 
small Kondo energy Tk ■ For calculations in the impor- 
tant temperature range T > Tk a solid treatment of the 
fermion constraint is therefore desirable. 

As far as the low-temperature behavior is concerned, 
the criticism of the auxiliary-particle approach often ex- 
pressed in view of the uncontrolled handling of the con- 
straint may be refuted on the basis of the results pre- 
sented here. However, one has to keep in mind that the 
above arguments are based on the assumption that the 
approximation method (whether based on self-consistent 
diagrams or functional integrals) does conserve the local 
fermion charge Qi . 

The Popov-Fedotov approach opens the way to us- 
ing resummed perturbation theory in specific strongly 
correlated systems, on the basis of standard Feynman 
diagrams, and for all temperatures. It requires identi- 
fying and performing the summation of physically rele- 
vant terms (diagram classes), which, however, remains a 
challenge for these systems. The self-consistent approx- 
imation presented here, for example, still fails to satisfy 
the notoriously hard to meet sum rule on the local spin 
moment. More elaborate resummation schemes are nec- 
essary to correct this and other deficiencies, the reward 
being a detailed description of the spin dynamics not ac- 
cessible by any other analytical method. 
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APPENDIX A: PROPERTIES OF GREEN'S 
FUNCTIONS IN THE POPOV-FEDOTOV 
TECHNIQUE 

In this appendix we consider thermal (Matsubara) 
Green's functions in the Popov-Fedotov scheme. For 
some operators A and B , which are both either fermionic 
(s = +1) or bosonic (s = —1), the Green's function is de- 
fined as^, 

G{t) = (~ s )(T t {A(t)B(0)})^ v , -p < r < [3 , 

with (3 — 1/T, the thermal expectation value as defined 
in Eq. (|2"U)) and (JT3J) , the Hamiltonian given by Eqs.©, 
([T|), and the usual "time" -ordering symbol 

T t {A(t)B(0)} = A(t)B(0)Q(t) - sB(0)A(T)e(-T) . 

The fact that H ppv , Eq.©, is non-Hermitian, does not 
influence the (anti-) symmetry properties resulting from 



the cyclic invariance of the trace. Therefore it is sufficient 
to consider r > , i.e., 

G(r) = { -s)^Trf[e-V-^Ae-^B] . 

Using Eqs. tfTTj) and © this becomes 

g m = (-)^ e e • 

c Q ,n Q c' Q ,n' Q 

■ (cq ,n Q \A\ c' Q , ri Q ) (c' Q , n' Q \ B \ c Q , n Q ) ■ 

e -PE{c Q ,n Q ) e T[E{c Q ,n Q )-E{c' Q ,n' Q ) ] _ 
N L N L 

.( J] e ifQ*) ( JJ e -*#[Q«-Q'.]) (Al) 

k=l q=l 

Q q and Q' q denote the auxiliary charge on lattice site q 
as it appears in the charge configurations cq and c'q , 
respectively. 

Physical propagator: The simplest physical Green's 
function is the dynamical spin susceptibility (fT5)) . (|19[) . 
for two lattice sites I and m , 

xZir) = {%{s^T)smr v ■ 

Due to the property (|13p of spin operators, the fermion 
charge on the sites I, m is automatically constrained to 
1 in Eq.jH]), i.e., Qi = Q m = Q[ = Q' m = 1 . For all 
other sites p ^ I, m , the orthonormal matrix elements 
in Eq. (|A1|) lead to Q p — Q' p , thus we have c'q = cq , 

and the second factor (fl^i e ) becomes 1. Now the 

term (Jl^i ex P(*f Qfc)) , in combination with the prop- 
erty (fl4|) of the energies, leads to a cancellation of all 
unphysical states with charge Q p — and Q p — 2 on any 
lattice site p ^ l,m . Thus, only physical states with a 
single fermion per site, cq = c'q = c P Q hys = (1,1, ...,1) 
remain in Eq. (|Aip . Using the notation 

E n = E{dQ h v\n Q ) , \n) = \J G h »%n Q ) 

for energies and states in the physical subspace, Eq. (|Al[) 
reads 

=^7E^i s ri^^i^i") e ^ £ " eT(E ' , " E '" ) 

n,n' 

(A2) 

With the result (fT?)) for the partition function, this is 
exactly the expression we would have obtained directly, 
working in the physical Hilbert space. 

Green's function of the fermions: The fermion prop- 
agator is not a meaningful physical quantity. How- 
ever, within a self-consistent diagrammatic expansion 
of, e.g., the dynamical spin susceptibility, the renormal- 
ized fermion Green's function is of technical importance. 
Therefore it is useful to derive some exact properties of 
the Green's function 

Gim.Ar) = ~(%{fUr)fLm) WV ■ (A3) 
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The matrix elements in Eq. (|Al[) now read 

(cQ^ol/i<r|cQ,nQ>(c / Q ,n / Q |/t l<7 |cQ,nQ) . 

f} na increases the auxiliary charge at lattice site m by 
1, which can only be compensated by fi„ , that is, the 
exact fermion propagator is local, I = m . For all sites 
p =/= I the arguments from above hold: The orthonormal 
wave functions lead to Q p = Q' p , and the unphysical 
constributions with Q p = 0, 2 cancel. The states that 
remain in the trace in Eq. (|Al[) then have Q p — 1 at all 
sites p I and some charge Qi,Q[ G {0, 1, 2} at site I . 
With the notation 

\Q,n) 1 = \c q ,n Q ) for Qi = Q , Q p = 1 , p ^ I , 

and similarly E l (Q,n) for the eigenenergies, the Green's 
function (|Aip . (|A3[) becomes 



G/ m>CT (r) = (-l)8h 



(0 



N L -1 



Q=Q,ln,ri 



(A4) 



■\XQ,n\f la \Q + l,n') 
•exp (Q,n)) exp (*|q) • 

•exp(r[£; / (Q,n)-£; i (g + l,n') 



'2/3 



]) 



The energies and states that occur in Eq. (|A4[) are now 
denoted by 

Q= 1 : E l (l,n)=E n , 

= 0,2: ^(0,n') = S'(2,n')=^ , 

and 

Q = l : |l,n)' = |n) , 

= : //j0,n')' = \a,n') 1 , 

= 2 : // CT |2,n') ! =±\-a,n') 1 . 

According to Eqs. (fTU)) . (fTT|) the £?„ and |n) are the 
eigenenergies and -states of the model Hamiltonian H , 
Eq.Q . Referring to Eas.ljTg]) and JSJ), the E l n , can be 
interpreted as the eigenenergies of H with a "defect" at 
site I, i.e., with all couplings J to the spin at site I set to 
zero. The states \cr,n'} 1 therefore contain the orientation 
a G {1,1} OI the resulting free spin at site I as a good 
quntum number. The set of quantum numbers n' as well 
as the E l n , do not depend on a . Note that the number 
of states is the same, #{|«)} = #{|a, n')} = 2 Nl . 
The fermion Green's function (|Aip now reads, 



GiAr) = 4£- 



(A5) 



• e-^"|H-a,n')f e T[ 



'\1\2t[E„-E 1 ,+i 



2r! - 



ie- fiE n'\(n\a, n') l \ 2 e 



In frequency space, the fermion propagator is given by 

G lm<a (iu) = / dre^ T G lm ^(T) , (A6) 
Jo 

with the fermionic (odd) Matsubare frequency u> = (2n + 
1)5. Inserting Eq. (|A5P into Eq. (|A6p and utilizing the 
earlier result Ea-lfTfj) for Z ppv , we obtain the fermion 
Green's function. 



Gi m . a (iu) = Si m / de - 



(A7) 



with the complex-valued spectral function 

G lt<T (£)=pUe)+iph(e) , (A8) 



Z 

\{n\a,n') l \ 2 8(e-[E n -E l n ,}) + 
+ \{n\-a,n') l \ 2 5(e-[E l n ,-E n }) 



n,n' 

.{\(n\-a,n')f5(e-[E l n ,-E n })- 
-\{n\o,n>) l \H(e-[E n -E l n ,])} 

Sum rule and symmetry: The \n) as well as the |cr, n') 
form a complete normalized basis in the physical Hilbcrt 
space, 

5>><n| = l , £V,"')M< = 1 > 

and therefore integrating Eqs. (|A8[ ) over e leads to the 
sum rule 

/CO />OC 
depL(e) = l , / d £ p^( E ) = , 



(A9) 



/ dEGi l0 .(e) = l. 



From Eq. (|A8p we can also read off a "particle-hole" 
symmetry, 

Pi,t( e )=Pi,i(- e ) ' ^t( £ ) = -^,i(- e ) . 



Gi, T (e) = G u (-e)* 



(A10) 
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In the paramagnetic phase, where the overlaps in Eq. (IA8 
are spin degenerate, 



"J 



(n| T,n') = (n| j,n') 



Eq. (|A10p simplyfies to the result already quoted in 
Eq.®. 

The expectation value (Qi) : In order to conclude this 
appendix, we discuss the average auxiliary charge (Qi) 
at site I . 

In the enlarged Hilbert space the expectation value 
(Qi) ppv can be formally calculated; however, although Qi 
is a gauge-invariant operator, it does not fulfill the prop- 
erty (fT5|) of physical observables, and therefore the result 
becomes meaningless. This is most easily demonstrated 
by explicitly calculating (Qi} ppv : 

Using the Green's function (|A3[) it may be written as 



(Qi 



ppv 



^G< )CT (t = 0_) 



In Eq. (|A5| the propagator has been given for r > , 
which can be utilized by help of the anti-symmetric prop- 
erty of fermionic Green's functions, 

r < : G(r) = -G(r + f3) . 

Thus we find from Eq. (|A5|) , setting r = (3 , 



-PE T . 



Z 

+ i 



(n\n) 



^J2J2e-^\a,n'\a,n'y 



n' cr 



l 



(Z + iZ l ) 



Here Z is the partition function of H , Eq. (fTJ) , in the 
physical subspace, while Z l is the partition function of 
H with the "defect" at site I , i.e., with all couplings J to 
the site I made zero. Since all interactions in H are short 
ranged, Z and Z l become equal in the thermodynamic 
limit. 



iV L ^oo: (Q,)PP" = (l + j_) _> (l + ») (All) 

Note that Z and Z' contain the same number of states, 
(2)^ . 



APPENDIX B: DERIVATION OF THE 
SELF-CONSISTENT EQUATIONS ([57]) 

In this appendix the intermediate steps in going from 
Eas.fiT]) to Eqs.gT)) are presented. 

Starting point is the spectral representation (|43|) or 
(|A"7|) of the fermion Green's function (|4"Te)) . The Mat- 
subara frequency ito can be analytically continued to the 
complex plane, ito — ► z , with G(z) showing a branch cut 



at Im(z) = — 57j • Close to this cut, at z — uj — =p i0+ , 
(0+ is a positive infinitesimal) we have 

G(u -t^T i0 + ) = ±i7rG(w) + G(w) , (Bl) 
2£ 

with the spectral function G and its Hilbert transform 

G(e) _ 

— e 



G(w) = P / de 



For the susceptibilities II(ii/) and D(iv) , appearing in 
Eq. (|41ap and ()41cj) , the usual analytic continuation of the 
bosonic Matsubara frequency iv to the real axis applies, 
iv — ► u; + z0+ , 



n(w) -> n(w + io+) = n'(w) + in"(w) , 

D(iz/) -> + i0+) = L»'(w) +iD"(w) . 



(B2) 



The imaginary part D"(uj) represents the spectral func- 
tion of the effective local interaction, 



1 f 00 D"(e) 
D(iv) = - de 

7T ./_„„ £ - iZ/ 



Note that D"(uS) obeys the symmetry 
D"(-u) = -£>» , 



(B3) 



(B4) 



which comes from %(q, tS) = x*(— q, — ui) in Eq. (l41cD . 
For Xl(iv) , equations similar to (|B3|) and (|B4p hold. 

The fermion self-energy (|41d[) is re-written using the 
spectral representations (|A7|) and (jB3|) . 



E(iw) = / de Z>"(e) / de'G(e') 



/? 4-^ (iv — e){iv + iuj + ifp — e') 
^- / de£>"(e) / de'G(e') 



+ + e — e' 



/ and <7 stand for the Fermi and Bose function. Appar- 
ently, E(iaj) obeys a spectral representation similar to 
Eq. (jA"7]) . namely, 



E(iw) = y de 



E(e) 



zcj + Ira — 6 



(B5) 



with the (complex valued) spectral function 

S(w) = ^ | d £j D"( e )G( e + c)[ ff ( £ ) + /( £ + c-z^)j 



and its Hilbert transform E , given in Eq. (|47ip above. 

Introducing the structure factor of the renormalized 
interaction, 

U(w) = [l+g{w)]D"(w) , 
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which is by Eq. (|B4p equivalent to 

U{-u) = g{uj)D"{uj) , 
and using the relation 

9 + f = g(i-f) + 0- + g)f 



the spectrum £ takes the form Eq. (|47h[) , with the short 
hands G+ and G~ defined in Eqs. (|47e|) and l[4"7fj) . 

The fermion spectrum G is obtained from the Dyson's 
equation (|41e[) using Eq. (|Bl[) . i.e., 



G(w) 



1 



1 



1 



w-E(w-^+iO + ) 
By inserting the decomposition 



S(a; — i— =F i0+) = ±i7rE(a>) + E(w) , 
Zp 

which results from the spectral representation (|B5[) and 
Eq. (|47i[) . we obtain G as given in Eq.( |47g[ ) . 

In the fermion bubble H(ii/j , Eq. (|41a() . the spectral 
representation (|A7|) of the fermion Green's function is 
inserted, and we arrive at 

Note that the imaginary-valued chemical potential //-^ = 
cancels in the denominator, since II represents an 
observable susceptibility. Apparently, U(ii/) obeys the 
usual spectral representation, similar to Eq. (|B3p . with 
the imaginary part 



H» = - / AsG{e)G{e-u) 



(B6) 



•[/(*-*- 

and the corresponding real part II' is computed via 
Eq.(l47cl) . 

It is convenient to introduce a structure factor for the 
bubble, 

s°H = [i + sH]n» , 



and with the relation 

[f(x-y)-f{x)][l + g(y)\ 



l-f( x )]f(x-y) , 



which is valid for arbitrary complex numbers x , y , we 
have 

s°M = | | dE^GMiwHl)]/^!) 



Using the notation G+ , G~ introduced Eqs.(|47e[). (|47fj) . 
the result stated in Eq. (|47d[) immediately follows. 

In order to compute II" (w) from S°(oj) , Eq. (|47b[) is 
used, which is a consequence of the symmetry II" {—to) = 

-n"(w). 

The last equation to be derived in this appendix is 
Eq. (|47a[) for the effective interaction U(u>) . Performing 
the analytic continuation iv — ► (u> + *0 + ) in Eqs. (|41c[) . 
(|41b[) . and using the decomposition (|B2[) . we find 



1 + J(q)n(w + iO+)| 



Eq. (|47aj) is now obtained using the definition of U (u>) and 
S°(u>) given in this appendix and the density-of-states 
jV(e) introduced in Eq. (|4"g]l . 

Particle-hole symmetry: In the Appendix |A] above, 
a symmetry for the spectrum G of the fermion Green's 
function has been derived in Eq JAlOj ). namely 

G(-uj) = G{uj)* . 

Accordingly, the spectra G + and G~ introduced in 
Eqs. (j47e|) . ijiTfj) obey the relation 



G-(-w) = G+(w)* 



(B7) 



This may be used to simplify the Eqs. (|47|) somewhat by 
eliminating G~ : Applying Eq.(|B7jl to Eq. (|47hj) leads to 



E(w) = — / def/(e)[G + (w"£) + G+f-w-e)*] . 

47T J 

(B8a) 



For S 10 , we start from Eq. (|47dj) by writing the expression 
twice and using the symmetry (|B7[) in the second term, 



S°(w) = -y d£[G+( £ )G-(£-w) + G-(-e)*G+(cj-e)*] 

By renamimg e — » (w— e) in the second term and applying 
Eq. (|B7|) once more, it follows 



S°(cj) = ~j deRc{G + (e)G + (uj-e)*} . (B8b) 

The remaining equations in (|47[) stay unchanged, except 
that Eq. (|47f|) becomes obsolete. 



APPENDIX C: THE SELF-CONSISTENT 
EQUATIONS USING REAL SPECTRA 

For a direct comparison of the self-consistent equations 
with those derived within average projection in Ref. 
(Eqs.(Al) in that reference), we find it instructive to re- 
write the Eqs. (|47p entirely in real- valued spectral func- 
tions. To that end we decompose all unphysical spectra 
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into real and imaginary parts as follows, 

G(u) = Pi (uj) + ip 2 (Lu) , (Cla) 
G+H = p+{uj)+ip+{uo) , (Clb) 

E(w) = &i(v)+ia 2 (u) , (Clc) 
E(w) = cti(w) +«ct 2 (w) . (Cld) 

Inserting these definitions in Eas.(|47p. with Eas. (|47hp 
and (|47d|) replaced by Eos. ([B8ap and (|B8b|) . we find the 
set of equations (|C2[) stated below. For completeness, in 
(|C2[) we also quote those equations from (|47)l . that re- 
main unchanged by using (|CJ1[) . Making use of Eq. ([3"7| 
the result reads 

7V(e) e 2 



f7(w) 


= 0V) / 


n"(w) 




n'(w) 






S°(w) 


-§/* 





|l + en(w + iO+)| 



n"( e ) 



;i-/(2w)]p!H + 

;i-/(2w)]^H - 



p 2 (w) 



2 cosh(/3w) 
2 cosh(/3w) 



Re 



Im 



(C2h) 

{A+A- - B+B-) - i(B + A- + A + B_) 

-I 

in J 

+ p+{-u)-e) 



2 (w 

The short hands A± , B± are defined as 
A± = U) — ai(u>) =fc 7ro"2(w) 

-B± = (72 (w) ± 7TCTi(u;) 



APPENDIX D: CALCULATION OF THE 
AUXILIARY-CHARGE RESPONSE ([69]) 



(C2a) 

(C2b) 
(C2c) 

(C2d) 

(C2e) 

(C2f) 
(C2g) 



deU(e)[p+(uj-e) + 
■f pf(-w-e)] 

A | d£ff( e )[p+( W - e )- 



(C2i) 
(C2j) 







p+(-w-e) 






o-i(e) 






u; — e 




'de 








uj — e 



(C2k) 
(C21) 



Eq.flTOD in the form 

r(q,iw) = 1 + ^2 [A{iu),iwx) 



(Dl) 



B(q; iw,icji)] T(q, iwi) 



with the matrices 



A(iu, iu>{) — —TG(iu>i) 2 D(iuj — iu>i) 

and B(q; ioj, iu>x) , the latter can be read off Eg. ([70)) . For 
the constituents of A and B we observe the symmetries 

G{-kS) = -G{iuS) , D(k,-iu) = D(k,iv) , 

which implies 

A{— iui, —iuii) = A(iu},iu>i) (D2a) 

and 



B(— iui,iu>i) = —B(iuj,iLUi) 
B(iLU,—iuJi) = —B(iuj,iLUi) 



(D2b) 



In this appendix the intermediate steps in going from 
Eg. ([70)1 to Eq. fTT)) are explained. We start by writing 



Here and in the following the wave vector q is not written. 
The vertex function is split into components behaving 
symmetrical (+) or anti-symmetrical (— ) under iui — > 

— ILU , 

corresponding to 

T( ± \iu J ) = ^[T(iu > )±T(-iuj)}. 

Using these definitions and the symmetry relations ()D2p 
in Eo. (|Dip . it can be seen that and decouple, 

T i+) (iLu) = l + ^A(^,^i)r (+) (iwi) , (D3a) 



(D3b) 

In the fermion bubble (|69[) only the symmetric compo- 
nent contributes (not writing q), 

x Q (o) = -2T^GM 2 r(+>M . 



Therefore, by identifying with V in Eg. ([D3ap . we 
arrive at Eg. (|7ip for the vertex function to be used in 
the fermion-charge response (|69p . 
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